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Abstract 

We consider spin-lattice relaxation processes for electrons trapped in lateral 
Si quantum dots in a [001] inversion layer. Such dots are characterized by 
strong confinement in the direction perpendicular to the surface and much 
weaker confinement in the lateral direction. The spin relaxation is assumed 
to be due to the modulation of electron ^-factor by the phonon-induced strain, 
as was shown previously for the shallow donors. The results clearly indicate 
that the specific valley structure of the ground electron state in Si quantum 
dots causes strong anisotropy for both the one-phonon and two-phonon spin 
relaxation rates. In addition, it gives rise to a partial suppression of the two- 
phonon relaxation in comparison to the spin relaxation of donor electrons. 



Typeset using REVTgX 



I. INTRODUCTION 



Recently, there is a growing interest in the physics of electron spin due to the enor- 
mous potential of spin-based devices. In these so-called "spintronic" devices, information 
is encoded in the spin state of individual electrons. Numerous concepts ranging from spin 
analogs of conventional electronic devices, to quantum computers (which utilize the Zeeman 
doublet of a confined electron clS 8b qubit)0 have been proposed. Electron spin states relax 
by scattering with imperfections or elementary excitations such as phonons. Hence, the spin 
relaxation time is a vital characteristic that determines the potential value of a spin-based 
device. 

In the present paper, we calculate the longitudinal spin-lattice relaxation time (Xi) for 
electrons confined in Si quantum dots (QDs). Being a relatively light semiconductor, Si is 
characterized by a weak spin-orbit interaction, which basically determines the strength of 
spin relaxation. In addition, only a small fraction of the isotopes in natural Si possess a non- 
zero nuclear magnetic moment. As a result, the electron-nuclear spin-flip process is expected 
to be slow. These two properties of Si, along with the tempting possibility of integrating 
the "quantum" part of a computer with the well-developed Si "classical" electronics, make 
Si an attractive material for spin devices. 

We concentrate on a particular design of QD based on a [001] inversion layer formed at 
the interface of Si and SiC>2 [see Fig. 1(a)]. The lateral confinement of electrons is assumed 
to be due to the attractive potential applied to the gate electrode deposited on top of 
the oxide layer. Alternative QD design based on Si/SiGe heterostructures is also possible. 
However, the SiGe-based design may face additional complications for quantum computing 
applications. For example, electron confinement in Si/SiGe heterostructures is relatively 
weak and the penetration of electron wavefunction to the SiGe barriers is inevitable. Since 
spin-orbit interaction in Ge is stronger than in Si, these structures are expected to have 
higher spin relaxation rates. 

In our calculations, we take advantage of the results obtained for spin relaxation of 
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shallow-donor electrons in SiBtofra In Refs. ^, |5], |6], the spin-lattice relaxation was assumed to 
be due to the modulation of electron g-factor by the phonon-induced strain. The parameters 
of the effective Hamiltonian were estimated from microscopic models. Later, they were 
also determined by measuring the g-factor under a static strain.! Taking into account the 
structure of the electron states in lateral QDs, we calculate the corresponding spin relaxation 
time for both single-phonon and two-phonon processes. Simple analytical expressions are 
obtained for the two limiting cases, namely, for T <C Ez and T ^> Ez, where T is the 
temperature in energy units and Ez is the Zeeman doublet splitting energy. The former 
case is more relevant to the conditions for quantum computation, while the latter is close to 
the typical conditions of EPR measurements. For the resonance frequency about 50 GHz 
and low enough temperature, the spin relaxation time is found to be several minutes with 
the single-phonon process providing the prevailing contribution to relaxation. For elevated 
temperatures, the two-phonon process becomes significant and the relaxation time can be 
substantially smaller. We also predict strong anisotropy in the relaxation rate. In particular, 
for the [001] inversion layer, the spin relaxation is suppressed if the applied magnetic field 
is parallel to the [001] or [110] directions. 

The rest of the paper is organized as follows: In Section II, we introduce the Hamiltonian 
for spin-lattice relaxation and provide expressions for the single-phonon and two-phonon 
relaxation rates. In Section III, the asymptotic dependences for different temperatures as 
well as numerical results are presented. Finally, Section IV is devoted to the discussion of 
the obtained results along with other potential mechanisms of spin-lattice relaxation. 



II. EFFECTIVE HAMILTONIAN FOR SPIN-LATTICE RELAXATION 

For crystals which possess inversion symmetry like Si, there is no analog of conventional 
deformation potential for the spin-flip process, which is often called Van Vleck cancellation.@ 
This is a direct result of the requirement for the Hamiltonian to be invariant under a gen- 
eralized inversion transformation C = JK, where J and K are spatial and time inversion 
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operators. In Refs. |], [|, || the following effective C-invariant Hamiltonian describing the 
modulation of electron ^-factor by strain has been proposed: 

H g = Aijf.iUijBf.ai, (I) 

where A^i are the coefficients, Uij is the strain tensor, B^ are the components of the 
magnetic field, and <j[ are the Pauli matrices; here and below we assume summation over 
the repeated indices. H g is written in the basis of the Bloch functions corresponding to the 
bottom of the conduction band and can be used for calculations of the spin transitions for the 
electron states described within the effective mass approximation. Nonzero coefficients A are 
determined by the symmetry of the crystal using the method of invariants.! In particular, 
for the A point of the Brillouin band in a diamond-like crystal that corresponds to the 
conduction band of Si, there are eight invariants and the Hamiltonian for a [001] valley can 
be written as:! 

H f m = {A x a z B z {u xx + u yy ) + A 2 a z B z u zz + (2) 
A 3 a z (B x u xz + B y u yz ) + A 4 (a x B x + o y B y )u zz + 
A 5 (a x B x + a y B y ){u xx + u yy ) + A 6 (a x B y + a y B x )u xy + 
A 7 (a x u xz + a y u yz )B z + A s (a x B x - o y B y ){u xx - u yy )) . 

Here, we introduced the factor /xb/2 for convenience, \x& being the Bohr magneton. The 
absolute values of the coefficients A can be determined in principle by using a many-band 
effective-mass expansion of the electron wavefunction in a uniform magnetic field, similar to 
that used in Ref. [|. In practice, this cannot be accomplished since the required momentum 
matrix elements are unknown. However, some qualitative considerations are possible. Since 
the terms of the expression for (/-factor contain the energy gaps between the coupled bands in 
the denominator, coupling of the closest bands is expected to be the strongest. For Si, there 
is a A' 2 band which is close to the Ai conduction band. These bands merge at the X point 
in the momentum space but are not coupled by either a spin-orbit interaction or momentum 
operators. Therefore, Ai <-> A' 2 coupling is not manifested in the effective mass or g-factor 
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of unstrained Si. However, these bands are coupled by the deformation potential. From 
the character tables of the A point and the corresponding invariants, it is easy to conclude 
that this coupling is realized by the deformation potential term proportional to u xy . Since 
H g ccontains only one invariant proportional to u xy , it can be concluded that the major 
contribution to the Hamiltonian is 

#[° 01 ] = X -Aii B {a x B y + a y B x )u xy} (3) 

where we drop the index of the coefficient Aq. This argument was initially used by Roth.! 

To proceed with the calculation of the longitudinal relaxation time T 1; we need to describe 
explicitly the system under consideration as well as the characteristic energy and length 
scales. As mentioned in the introduction, we consider Si lateral QDs formed at the Si/SiC>2 
interface, where the lateral confinement is due to the gate electrodes. For such a system, 
the lateral dimension of the QD ai at is typically on the order of one hundred nm, exceeding 
considerably the inversion layer thickness a 2 D- The energy structure of the electron levels is 
determined by the following parameters: (a) quantization energy in the 2D inversion channel 
E 2 d ~ ^ 2 /(2m2Dfl2D); 3 ) lateral quantization energy Ei at ~ h 2 /(2mi at a 2 at ); (c) intervalley 
splitting energy between [001] and [001] states, A; and (d) Zeeman energy E z = HbqB, 
where g is the effective g-factor of the confined electrons depending, in general, on the 
direction of the magnetic field. Here m 2 £> and m/ at are the electron effective masses in the 
direction normal and parallel to the inversion layer, respectively. For the [001] inversion 
layer, they just correspond to the longitudinal and transverse effective masses of Si. In the 
following, we assume that the conditions E z <C A, E [at and E[ at <C E 2 d are satisfied. In 
fact, this is a necessary requirement for spin qubit operation. Using modern technology, 
Ei at can be made about a meV or even higher. In contrast, A can be controlled in a 
much lesser degree, and is roughly proportional to the confining electric field in an inversion 



layer. Experimentally, A was measured to be up to 1 meV in strong fields (see Ref. |f0| ). 
In Fig. 1(b), we show schematically the energy levels and the electron transitions under 
consideration. The numbers "0" and "1" mark the levels of lateral electron confinement. 



The signs + and — denote the valley-split electron states, and "up" and "down" are the spin 
states. We assume that only the lowest Zeeman doublet can be populated, which means 
that the temperature T is much less than A and E\ at . The longitudinal relaxation time T\ 
is determined as Tf = W up -down + Wdown-u P , where W are the probabilities of spin-up to 
spin-down and spin-down to spin-up transitions. The solid arrows indicate the single-phonon 
transitions, while the dashed arrows correspond to the two-phonon process. Later in this 
section we comment on the possible two-phonon processes shown in the figure. 

In Refs. |] and ||], an alternative Hamiltonian to Eq. (H) was considered. It originates 
from the coupling of the donor singlet and doublet states by an applied magnetic field. It 
arises due to the different valley structure of these states and the anisotropy of the g-factor 
in the individual valleys. Although this Hamiltonian does not involve A <-> A' 2 coupling, 
its contribution is high because of a very small gap between the singlet and doublet states. 
However, there is no such mechanism in the [001] lateral QDs, which is probably the main 
difference with the case of spin relaxation for donor electrons. The reason for that can be 
easily seen. For the [001] inversion layer, the ground electron state is a combination of [001] 
and [001] states: 



where we dropped the spin index of wavef unctions. Here \ is the envelope wavefunction, -0ooi 
and V'ooi are t ne Bloch functions corresponding to the [001] and [001] valleys, respectively, 
and F are the coefficients which determine the valley splitting. Particular expressions for 
F can be found using a microscopic model, for example, that of Sham and Nakayama.0 
In our case we do not need explicit expressions for C. It is enough to use the fact that 
Zeeman Hamiltonian has identical forms for the [001] and [001] valleys and its intervalley 
matrix elements are zero. Therefore, the matrix element between and which is 
proportional to the overlap between them, is zero because these states are orthogonal. 

The Hamiltonian of Eq. ([J) is written in the representation where the basis electron 
wavefunctions correspond to the definite spin projections on the z axis. For calculation 




(4) 
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of Ti, it must be rewritten by using a representation with the definite spin projection on 
the direction of the magnetic field.0 This can be done following the standard procedure of 
transformation for Pauli matrixes under rotation.!^ Finally, we obtain the expression for the 
one-phonon relaxation rate, 1/Tj as: 
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4^ = !L^_^_(i + 2 iV r )sin 2 ^cos 2 20 + cos 2 ^sin 2 20)V / dttf [ x y \ y x } . (5) 
T W 4 g 2 p Y J s i 

In this equation, g is the slightly anisotropic g factor of confined electrons, / = gpsB /(2itK) 
is the resonance frequency, 9 and are the spherical angles of the magnetic field, p is the 
material density, Nt is the Planck phonon population for the energies equal to Ez, Q q is a 
solid angle in the phonon wavevector space q, e and n are the phonon polarization vector 
and the unit vector parallel to q, respectively, s is the sound velocity, and the summation 
is over the acoustic phonon branches. Equation (|5]) assumes that the phonon wavelength 
corresponding to the energy Ez is much greater than the lateral dimensions of QD. In this 
case, the form factor of the electron-phonon interaction is equal to unity and r/ 1 does not 
depend on the particular shape of the lateral confining potential. To check the validity of 
this approach, we performed calculations of T± assuming parabolic lateral confinement and 
found that for the lateral level separation of 1 meV, the obtained correction is less than 10 % 
even for / = 50 GHz. An additional assumption of the bulk-like phonon spectrum is made 
for simplicity. This probably introduces a greater error since the lateral QDs are normally 



situated close to the surface. As shown in Refs. |TJ, [15|, [16], the phonon modes in this case 
are essentially rebuilt, due to the interference of the incident and reflected phonons as well 
as the origination of Rayleigh waves, which strongly modify electron-phonon coupling. 

Equation (|5|) predicts strong anisotropy of the relaxation rate. If the magnetic field 
is parallel to the [001] or [110] directions, goes to infinity. This is because we used a 
single-parameter Hamiltonian of Eq. (|3|) . When rewritten for the basis with the definite spin 
projection on the direction of the magnetic field, it vanishes for these particular orientations. 
Of course, in experiments the relaxation for these cases is not expected to be suppressed 
completely, since the remaining terms of Eq. (||D have non-zero contributions to the relaxation 



rate. However, a significant decrease is expected. 



Assuming an isotropic acoustic phonon spectrum, we obtain 
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(1 + 2N T ) sin 2 (cos 2 20 + cos #sin 2(f)). 
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Here we take into account only TA phonons, which provide the major contribution to the 
relaxation rate. 

Let us now turn to the calculation of the relaxation time due to the two-phonon tran- 
sitions . This is a second-order transition where the electron is virtually scattered first 
to an intermediate state and then to the final state. One of the virtual transitions is ac- 
companied by spin flip, while the other occurs with spin conservation. The probability of 
a two-phonon transition can be found using the second-order perturbation theory (see, for 
example, Ref. [17]). For the spin relaxation of donor electrons, the intermediate electron 
state is represented by the excited doublet. In contrast, for the case of [001] lateral QD, 
the valley-split state cannot serve as an intermediate state. This is because the valley-split 
states are not coupled by the Hamiltonian of Eq. (Q). This can be shown using a argument 
similar to that applied for the proof of the absence of g-factor modulation due to coupling of 
the valley-split states. Other possible transitions are through the excited states of the lateral 
confinement. Since the intervalley splitting is controlled by the electron confinement in the z 
direction rather than in the lateral direction, these transitions are actually possible between 
the states having the same valley structure [see Fig. 1(b)]. If the phonon wavevector is very 
small, then the probability of such a transition goes to zero because the overlap functions of 
"0" and "1" states are orthogonal. To determine the probability of two-phonon transition 
for a finite phonon wavevector, we need to know the overlap function x explicitly. For cal- 
culations, we assume parabolic lateral confinement. In this case, the "lateral" electron state 
is determined by the two quantum numbers, l x and l y . The "0" state corresponds to l x = 0, 
l y = 0, and there are two degenerate "1" states with l x = 0, l y = 1, and l x — 1, l y — 0. 
Using the expressions for the wavefunctions of harmonic oscillator ,0 it is easy to obtain the 
necessary form-factor J: 
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J S I Ax dyd> XmXm exp + „ + ,.,)) = -L| exp , (7) 

where the subscript of x represent the l x and l y quantum numbers and k is expressed 
through the energy gap 8 between "0" and "1" states and the lateral effective mass: k = 
V m iat8/h. Here we take into account that the thickness of the inversion layer is much 
less than the typical phonon wavelength. In the following, we also assume that the typical 
phonon wavevector is less than k, and drop the exponent in the expression for J. Note, 
that J can be modified due to the diamagnetic influence of the magnetic field as well. In 
particular, it can lift the degeneracy of "1" states. In our calculations, we do not consider 
this effect. 

We can distinguish three contributions to the two-phonon relaxation rate: (a) due to 
emission of two phonons, (b) due to absorption of two phonons, and (c) due to phonon 
scattering. For each of them the rate can be expressed in a uniform manner: 

if) = 34 ' 10 ( m) 2 w£L sm2 9(cos2 24> + cos2 6 2m - (8) 

Here only the transverse phonons are taken into account, which provide the major contri- 
bution to the relaxation rate. We also assume that typical phonon energies are considerably 
less than 5. The spin-conserving virtual transition is treated within the deformation model, 
which for the [001] valley provides the interaction E = E^Ua + E 2 u ZZ) where E x and E 2 are 
the deformation potential constants. The coefficients D depend on temperature and they 
are different for each of the three processes mentioned: 

D em = /; dx *•(! - *)• (l + eMx ) t) _ J (l + exp((1 _ 1 , )/i) _ 1 ) , (9) 

D ab = / dxx 5 (l - xf — — — , 

Jo y ' exp(x/t) - 1 exp((l - x)/t) - 1' 

20 / poo if i \ 

D - = 17 U "* Al + Xf ^/t) - 1 V + exp((l + ,)A)-l j + 

r-oo j 
/ dxx 5 (x+ l) 5 r-rr 

Jo V ; exp((x + l)/t) - 1 



1 ' 



exp(x/t) — 1/ / ' 

where t = T/Ez- One can see, that the two-phonon relaxation rate is characterized by the 
same anisotropy as in the one-phonon relaxation rate. In the following section, we analyze 
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the asymptotic dependences of the relaxation rates for different temperature regimes and 
perform numerical calculations. 

III. RELAXATION RATES FOR LOW AND HIGH TEMPERATURES: 

NUMERICAL RESULTS 

For the limiting cases of T Ez and T 3> Ez, the relaxation rates obey simple power 
laws as a function of resonance frequency and temperature. In particular, ~ /™ 5 in the 
former case and 7^ ~ j-4y-i j n latter case. This is similar to the case of donor spin 
relaxation, ill 

The two-phonon relaxation time for T -C Ez is proportional to / 13 with the 
main contribution from the two-phonon emission. For T ^> Ez, the relaxation is mainly 
due to the phonon scattering and ~ /~ 2 T~ n . For donor electrons, Roth obtained a 
different power law.i This is because for donor electrons the form-factor of singlet-doublet 
transition in the lowest approximation is equal to unity, in contrast to Eq. ([?[) for the lateral 
QD where the form-factor is suppressed for long-wavelength phonons. 

In Figs. 2 and 3, we show the results of numerical calculations for one-phonon and two- 
phonon rates. We assume a magnetic field parallel to the [100] direction, p = 2329 kg/m 3 , 
the transverse sound velocity s t = 5420 m/s, E 2 = 10 eV, rrii at = 0.19 m , 5 = 2 meV , and 
g = 2. The coefficient A can be determined by the measurement of the donor g-factor in 
strained Si. With this method, A = 1.32 was obtained in Ref. 0. The dependence of the spin 
relaxation rate on the resonance frequency is presented in Fig. 2 for several temperatures. 
One can see that for the transition between the described power dependences at the 

(2) 

different temperature regimes is quite fast. This is not the case for T x , which is due to the 
different temperature dependence of phonon absorption, emission, and scattering rates. 

In Fig. 3, we plot the relaxation rates as a function of / under the condition T/Ez = 
const. This is relevant to the case of QD-based qubit, where this ratio must be kept small 
to ensure initial state preparation of the qubits. We see that under this condition the major 
contribution is provided by one-phonon scattering. 
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Of course, the predicted for small resonance frequencies huge relaxation times hardly can 
be measured in experiments. This is similar to the case of donor spin relaxation,! where 
actual experiments were undertaken for / about tens of GHz. 

IV. DISCUSSION 

Let us first summarize the major distinction of the spin-relaxation process for lateral QDs 
and shallow donors. First, we predict that the relaxation for lateral QDs is more anisotropic 
than that for donor electrons. In particular, for [001] inversion layer the relaxation rate is 
suppressed for a magnetic field parallel to the [001] or [110] direction. Second, the two- 
phonon relaxation in lateral QDs is, in general, weaker than for donor electrons and is 
characterized by different power dependences on the resonance frequency and temperature. 
Both of these features arise due to the different valley structure of the electron states in 
lateral QDs in comparison to that of donor electrons. 

We have to stress that the first conclusion relies on the model used for calculations, which 
must be checked by experiments. This is because the values of the coefficients A n in Eq. (0) 
are determined not only by the energy gaps between the bands, but also by a number of 
interband matrix elements, which are unknown. Strictly speaking, the experiments with 
donors can not be considered as a rigorous proof of single- valley Hamiltonian of Eq. (|j) . In 
fact, for the donor electrons, eight invariants of Eq. (§) are transformed to three invariants 
after summation over the valleys. The term measured in Ref. §, H ~ u xy (a x B y + a y B x )) + cp 
where cp stands for cyclic permutations, is obtained from several terms of Eq. @, not only 
from that of Eq. @ . 

There is another mechanism that can modify the spin relaxation in the lateral QDs. This 
can take place if some quantum levels originated from the longitudinal [001], [001] valleys 
and the transverse [100], [100], [010], [010] valleys come close. In this case, the intervalley 
coupling can mix these two groups. Mixing of these states and the ground state of the QD by 
both steady-state and strain-induced contributions of Zeeman Hamiltonian can be possible. 



11 



As a result, both one-phonon and two-phonon relaxation will be modified. According to self- 
consistent calculations,^ such a situation is possible for particular parameters of Si inversion 
layer. 

Finally, it is necessary to take into account alternative forms of the Hamiltonian which 
also cause spin relaxation. For example, the following Hamiltonian is possible: 

H sp = Dijki-^—cTi. (10) 

As with H g , this Hamiltonian is invariant under inversion C . Physically, it describes the 
splitting of the spin states by nonuniform deformation. A similar relaxation process pro- 
portional to the third power of the phonon wavevector, was discussed previously in Ref. 
|18| . We believe that the disagreement between the theoretical calculations of 7\ based on 
H g and experimental measurements! can be partially explained by the contribution of H sp . 
Indeed, it is easy to check that H sp provides the same dependence of the one-phonon rate 
on frequency and temperature as H g . This is because the energy conservation makes the 
magnetic field and the phonon wavevector proportional to each other. In Ref. 0, the mea- 
sured one-phonon rate was about two times higher than the calculated one. The two-phonon 
relaxation was found to obey the predicted temperature dependence, but in contradiction 
to H g , demonstrated no dependence on the magnetic field. Assuming that the contribution 
of H sp to one-phonon relaxation is about that of H g , the contribution of H sp to two-phonon 
rate in this case can be much greater than that of H g since the typical phonon energy for 
two-phonon process is much greater than that of Zeeman splitting. Hence, this assumption 
can explain both the quantitative disagreement obtained for one-phonon relaxation and the 
qualitative one obtained for two-phonon relaxation. Of course, this idea must be checked 
experimentally. In particular, important information can be obtained from the measure- 
ment of possible spin splitting in nonuniformly strained Si. Consideration of the symmetry 
properties and the absolute value of such splitting will be provided elsewhere. 

Finally, we would like to stress that the electron spin relaxation in QDs is of the same 
order of magnitude as that obtained for donors. It is much longer than the T\ in III-V 
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compounds, which proves a good perspective of Si for quantum information processing. 
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FIGURES 

FIG. 1. (a) Schematics of the lateral QD. Confinement in the [001] direction is achieved as 
in the conventional inversion layers, while the lateral confinement is provided by the electrodes 
(black boxes), (b) Schematic illustration of energy levels in lateral QDs. (0) and (1) mark the 
different states of lateral confinement. Valley-split + and — states are further split by the magnetic 
field into spin-up and spin-down states. The solid (dashed) arrows show electron transitions under 
one-phonon (two-phonon) relaxation. 

FIG. 2. Dependence of one-phonon (solid lines) and two-phonon (dashed lines) spin relaxation 
rates on resonance frequency for temperatures of 0.5 K, 1 K, and 2 K. 

FIG. 3. Dependence of the spin relaxation rates on resonance frequency for the case where 
the ratio of temperature and Zeeman splitting t is kept constant. The one-phonon rates for the 
considered values of this ratio cannot be resolved on this scale and is presented by one curve. 
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